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Abstract 


We investigate the chemical equilibration of the parton distributions in col¬ 
lisions of two heavy nuclei, assuming the partonic fluid to be ideal as well 
as viscous. The initial conditions are taken from HIJING calculations for 
Au-|-Au collisions at RHIC and LHC energies. It was seen that when the vis¬ 
cous drag is taken into account in the fluid flow, the life time of the plasma is 
increased by nearly a factor of 2. The temperature as well as fugacities evolve 
slowly than their ideal counterpart. The photon and lepton pair production 
was also investigated. There is a two fold increase in the photon and lepton 
pair numbers with viscosity on. The increase in the large pT photons and 
the large invariant mass lepton pairs are due to slower rate of temperature 
evolution. 
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I. INTRODUCTION 


In relativistic heavy ion collisions, a strongly interacting matter called quark-gluon 
plasma (QGP) is expected to be formed. The two colliding nuclei can be visualized as 
two clouds of valence and sea partons, which pass through each other and interact [Q]. The 
multiple parton collisions can then produce a dense plasma of quarks and gluons. After 
its formation, the plasma will expand, cool and become more dilute. If quantum chromo¬ 
dynamics admits hrst order deconhnement or chiral phase transition, it is likely that the 
system will pass through a mixed phase of quarks, gluons and hadrons, before the hadrons 
lose thermal contact at sufficient dilution and stream freely towards the detectors. 

Several questions about the structure of the matter formed in these nuclear collisions 
arise. Does the initial partonic system attain kinetic equilibrium? Probably yes, as the 
initial parton density is large and the partons suffer many collisions in a very short time 
1^. Does it attain chemical equilibrium? This will depend on the time available |^-[^ to 
the partonic system before it convert into a mixed phase or before perturbative QCD is no 
longer applicable. The time available for equilibration is perhaps too short (3-5 fm/c) at the 
energies (a/s < 100 GeV) accessible to the Relativistic Heavy Ion Gollider (RHIG). At the 
energies (y^ < 3 TeV/nucleon) that will be achieved at GERN Large Hadron Gollider(LHG), 
this time could be large (more than 10 fm/c). If one consider only a longitudinal expansion 
of the system, the QGP formed at LHG energies could approach chemical equilibrium very 
closely, due to higher initial temperature predicted to be attained there. The question of 
chemical equilibrium has the added importance as the preequilibrium phase may influence 
the yield of certain quark-gluon plasma signals, such as lepton pairs, photons and hadron 
containing heavy quarks. 

Space-time evolution of quark and gluon distributions has been investigated in the frame¬ 
work of parton cascade model [||], which is based on the concept of inside-outside cascade 
M and evolve parton distributions by Monte-Garlo simulation of a relativistic transport 
equation involving lowest order perturbative QGD scattering and parton fragmentations. 
Early calculations were done by assuming fixed pt and virtuality cutoffs for the partonic 
interactions to ensure the applicability of the perturbative expansion of QGD scattering 
process. From numerical studies P,pTHT5[| three distinct phases of parton evolution can be 
distinguished: (1) Gluon thermalise very rapidly, reaching approximately isotropic momen¬ 
tum space distribution after a time of the order of 0.3 fm/c. (2) Full equilibration of gluon 
phase space density takes considerably longer. (3) The evolution of quark distributions lags 
behind that of the gluons, the relevant QGD cross sections being suppressed by a factor of 
2-3. In the recently formulated self-screened parton cascade model a early hard scatter¬ 
ing produces a medium which screens the longer ranged color fields associated with softer 
interactions. When two heavy nuclei collide at sufficiently high energy, the screening occurs 
at a length scale where perturbative QGD still applies. This approach yields predictions 
for initial conditions of the forming QGP, without the need for any ad-hoc momentum and 
virtuality cutoff. This calculations also indicate that the QGP likely to be formed in such 
collisions are far from chemical equilibrium. 

Essentially same picture emerges from the study of Biro et al |^ . Their calculations were 
motivated to obtain a lucid understanding of the dependence of different time scales on var¬ 
ious parameters and model assumptions and a better physical understanding of the infrared 
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cut off required in the parton cascade model .The parton cascade model is a complex model 
and it is difficult to unravel these dependence from model calculations. They assumed that 
after a time Tjso, the evolution of partonic system is governed by the energy-momentum 
conservation law. They derived a set of rate equations describing the chemical equilibration 
of gluons and quarks, including the medium effect on relevant transport coefficients. The 
initial conditions of the partonic system was obtained from HIJING calculations. Their re¬ 
sults support the picture emerging from parton cascade model, that the plasma is essentially 
a gluon plasma. The gluon distributions attain near chemical equilibrium at LHC energies, 
and the quark distributions lag behind. At RHIC energies, the gluon as well as quark dis¬ 
tributions are far from equilibrium. In their calculation the partonic fluid is assumed to 
undergo longitudinal expansion only. Effect of transverse expansion on chemical equilibra¬ 
tion was studied in ref.[ [^,0]. At RHIC energies, results for pure longitudinal expansion 
and transverse expansion are quite similar except at large values of transverse radius. How¬ 
ever, at LHC energies, the transverse expansion changes the scenario, the plasma initially 
approach the chemical equilibrium, but then is driven away from it, when the transverse 
velocity gradients develop. Interestingly, total parton multiplicity as well as lepton pair and 
photon productions are not much affected by the transverse expansion. 

In all these calculations P, P!7| , P!B| , the hydrodynamic evolution of the partonic system 
was assumed to be ideal. It is well known that dissipative effects like viscosity can affect the 
flow . The system will cool slowly as viscous drag will oppose the expansion. As the 

chemical equilibration rate depend on the cooling rate, it will also change. In the present 
paper, we will explore the effect of viscosity on the evolution and equilibration of parton 
distributions. We will also investigate its effect on the lepton pair and photon production, 
typical QGP signals. 

The paper is organised as follows: in section 2, we set up the basic equations for hy¬ 
drodynamic and chemical evolution of the plasma, also define the viscosity coefficients. In 
section 3, the results will be discussed. In section 4, we study the effect of viscosity on 
photon and dilepton productions. Summary will be given in section 5. 


II. HYDRODYNAMIC EXPANSION AND CHEMICAL EQUILIBRATION 

A. Basic Equations 

We assume that after the collision, a symmetric partonic system (QGP) is formed in 
the central rapidity region. We also assume that the partonic system achieves a kinetic 
equilibrium by time Tiso, when the momenta of partons became locally isotropic [Q]. Lo¬ 
cal isotropy in momentum can be said to be established when variance in the longitudinal 
momentum distribution equals the variance in the transverse momentum distribution. This 
occurs after a proper time of about 0.7A/ |^, Xf being the partonic mean free path. At 
collider energies, this time corresponds to Tiso ~ 0.2 — 0.3 fm/c 0. Beyond further ex¬ 
pansion of the partonic system can be described by hydrodynamic equations. The approach 
to chemical equilibration is then governed by a set of master equation which are driven by 
the two body reactions( 5 f 5 f <->■ qq) and gluon multiplication and its inverse process, gluon 
fusion {gg <->■ ggg). The hot matter continue to expand and cool due to expansion and 
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chemical equilibration, till the temperature falls below the critical value (Tc = 160 MeV), 
when we terminate the evolution. 

We assume that the hydrodynamic expansion is purely longitudinal. As indicated in 

3,0 


the analysis 


at RHIC energies, the transverse expansion effect is minimal. It does 


affect the parton equilibration rate at LHC energies, the effect showing sensitive dependence 
on the initial condition of the plasma. For the initial condition as obtained from HIJING 
calculations, this effect is not large Ezia- 

For boost-invariant longitudinal flow, the energy-momentum conservation equation for 
the partonic fluid including the viscosity, is the well known equation EMU. 


de 

dr 


e + P — 


4/3rj+C 


= 0 


r 


( 1 ) 


where rj and ( are the shear and bulk viscosity coefficients. In the QGP, the bulk viscosity 
coefficient is zero. It is apparent that in one dimensional flow, the effect of viscosity is to 
reduce the pressure. One must note that, under no circumstances, the viscous drag can 
exceed the pressure. 

Equation of state need to be specihed for solving the hydrodynamic equation |l|. For 
partially equilibrated plasma of massless particles, the equation of state can be written as 

i 


e — 3P — [a 2 Xg -l- b 2 {Xq + Xg)]T^ (2) 

which imply a speed of sound = 1/a/3. In Eq.||, 02 = Svr^/lh, 62 = Tvr^A^j/dO, Nf ^ 2.5 
is the number of dynamical quark flavor. A* is the fugacity for the parton species i. Here 
we have defined the fugacities through the relations: 


'^g 

Tlq = Xqflq 
TLq = Xqflq 


(3a) 

(3b) 

(3c) 


where hj is the equilibrium density for the parton species i: 

16 

hg= ^C(3)T3 ^ 

^C(3)/V/T» = (4b) 


We further assume that Xq = Xq. Denoting the derivative with respect to r by an overdot, 
eq.^ can be rewritten as, 

^ _ I _= 0 15') 

Xg -|- b(Xq -|- Xq) T 3r 3 \fi 2 Xg -|- & 2 (Ag -|- Aq)]T^ 

where 5 = 62/^2 = 21A^j/64. If the partonic fluid is considered to be ideal fluid, 77 = 0 and 
the above equation can be integrated to obtain. 
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[Ag + b{Xq + Xq)]^^‘^T^T = const. (6) 

For a fully equilibrated plasma, (A^ = A^ = Ag = 1), and we obtain the Bjorken scaling 
law T^t = const. For viscous fluid, eq.^ can not be integrated and has to be evaluated 
numerically. 

The master equation for the dominant reaction channels gg ^ ggg and gg ^ qq for 
chemical equilibration are. 


>3 .^3—> 2 ) ^g^g^g') i'^) 

dfj.i^nqU^') d^i^nqU^') {jlgRg^q nqRq^g') ( 8 ) 

in an obvious notation. Using Juttner distribution for the phase space distribution for the 
partons, the rates can be factorized as |Q], 


ng{R2^3 - i?3^2)= ^o-3’^g(l 
1 2 / 

^g^g^g ~ ^qRq^g= 



Uqnqn'^g 

hgfiqnl 


) 


(9a) 

(9b) 


where and (T 2 are thermally averaged, velocity weighted cross sections, 

c^3=< (T{gg ggg)v > (lOa) 

cr 2 =< a{gg qq)v >, (10b) 


Dehning density weighted reaction rates; R 3 = \cr 3 ng and R 2 = \o' 2 ng, eq.|| can be 
written as. 


^ + 3— + i?3(l - Xg) - 2i?2(l - 


X 


T 


Ul /Xq Xq . 


' +3T + -=^2r^)r -A) 

T T bi Xq Xq 


( 11 ) 

( 12 ) 


Biro et al p] have evaluated the reaction rates R 2 and R 3 , taking into account the color 
Debye screening and the Landau-Pomeranchuk-Migdal effect suppressing the induced gluon 
radiation. The results are. 


/? 2 ~ 0.24A^ja^AgT ln(1.65/asAg) 

R 3 = 2 .W,T{ 2 Xg-Xiy/^ 


(13) 

(14) 


The coupled Eqs(^^|^ determine the evolution of T(r), Ag(r) and Aq(r) towards the 
chemical equilibrium and are solved numerically. It may be noted that the rate equations 
(^,^) may give negative growth rate for the fugacities, in case of very slow cooling. For 
viscous fluid, viscosity being temperature dependent, the cooling is slowed down, and such 
a situation may arise at early times when the temperature of the partonic system is very 
large. However, physics of the problem require that fugacities can grow only. Thus, while 
we impose the additional constraint that. 


solving eqs(pT 


Xg^g > 0 


(15) 
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B. Viscosity coefficient 


Several results for viscosity coefficients for QGP system is available |ll^J^,|22|-p4| . Thoma 


241 has calculated the shear viscosity coefficient of the QGP, in the relaxation time approx¬ 


imation. Screening effects were taken into account by using an effective perturbation theory 
developed for the finite temperature QGD, in the weak coupling limit |^. His results agree 
well with the variational calculation of Baym et al |2^ 
calculation. 


and will be used in the present 


Vq 


Vg 


T3 

0.82——-—— 

ln(l/as) 

T3 

0 . 20 -—-—— 

ln(l/as) 


(16a) 

(16b) 


For plasma, away from equilibrium, the shear viscosity coefficient for QGP then can be 
written as. 


I) = + A,)), 


(17) 


with r]g^q given in eq.|TC. 


We have assumed that the dominant reactions towards the chemical equilibrium are the 
99 ^ 999 and gg ^ qq reactions. The reaction rates for these processes (eq.|^|T^ can be 
used in the kinetic approximation to obtain the shear viscosity coefficients. In the kinetic 
approximation. 


V ~ —rii <pi> \i 

15 


(18) 


where rii and < pi > are the density and average momentum of the particles of type i and 
Xi is its mean free path. Using < p >= 3.2T the viscosity coefficient can be written as. 


12.8 aiT^ 

30 + 


(19) 


In the following calculations we will use both the viscosity coefficients eq. o and eq.|T9| 
and will be referred as vicosity coefficient I and II respectively. 


C. Results for RHIC and LHC energies 

The initial conditions for the hydrodynamic evolution are listed in table (1) P]. They 
are the results of HIJING model calculation, which is a QGD motivated phenomenological 
model, as only initial direct parton scatterings are taken into account. Thus there are some 
uncertainties in these parameters. However, they suffice our purpose of demonstrating the 
effect of viscosity on parton equilibration process. 

In fig.la,b and c, we have shown the evolution of the temperature, the gluon fugacity 
and the quark fugacity at RHIG energy. The solid line corresponds to the ideal flow, i.e. no 
viscous drag, the dashed and the long dashed lines are with viscous drag, corresponding to 
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shear viscosity coefficients I and II As expected, with viscosity, temperature of the partonic 
system evolve more slowly. Viscosity oppose the expansion, consequently, cooling is slowed 
down. The lifetime of the QGP phase is considerably increased in presence of viscosity (by 
nearly a factor of 2). It is apparent that this will have considerable effect on signals of QGP 
such as dileptons and photons. We also note that temperature evolution do not differ much 
with different prescriptions for r]. Life time being changed by less than 2%. 

The evolution of gluon fugacity is much more interesting. With viscosity, gluon fugacity 
also evolves more slowly. Viscosity impedes the chemical equilibration process. We hnd that 
upto 2 fm, Xg donot change, it then increases slowly, much more slowly than for the ideal 
fluid. This is contrary to the expectation from ideal fluid flow (eq.^ that if the temperature 
evolve slowly, fugacities will increase at a faster rate. Indeed for viscous fluid flow, eq.|| is 
no longer valid. The behavior of Xg can be understood by recasting eq.|| as. 


W_3^_V«3(1 


Xg) — 2i?2(l 




( 20 ) 


Then as the temperature decreases more slowly, rhs of eq.^ will decrease as a result 
of which, fugacity will also evolve slowly. Thus viscosity, affect both the temperature and 
fugacity evolution, requiring them to evolve more slowly than their ideal counterpart. We 
also note that in the viscous partonic system, the gluon fugacity could not attain the value 
reached in the ideal case, by the time the system is cooled to Tc, even though the life time 
is considerably increased. 

In hg.lc, the quark fugacity is shown. They show a similar behavior as gluon fugacity. 
As with gluon fugacities, quark fugacities also evolve much more slowly if the viscosity is 
turned on, for the same reasons as described above. By the time is reached, they are far 
behind chemical equilibrium. Here also, the fugacities do not attain the value reached in 
the ideal case at Tc. 

The results for evolution of temperature and fugacities at LHG energies are nearly the 
same, as shown in fig.2a,b and c. Here also, with viscosity, plasma cools slowly than the ideal 
fluid, life time of the plasma increasing by a factor of 2. The gluon and quark fugacities also 
evolve slowly. In contrast to RHIG energy, gluon and quark fugacities attain larger values 
at Tc. This is due to their comparatively large initial values. Thus our initial conjecture 
that chemical equilibrium is attained at LHG energies but not at RHIG is substantiated. 
However, we note that even at LHG energies, the plasma is not fully equilibrated, Xg ~ 0.7 
and Xq ~ .5, and are quite far from equilibrium. 


III. THERMAL PHOTONS AND LEPTON PAIRS 
A. Photon spectra 

Thermal photons and lepton pairs are primary probes of quark-gluon plasma. Being 
weakly interacting, they carry the information of the early phase of the hot fireball created 
in the collision. Thermal photons from QGP has their origin in the Gompton {qg qj) 
and annihilation processes {qq —> gj) processes. For plasma away from equilibrium, their 
rates have been calculated 
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( 21 ) 


^ dN 2 aas, ,^99 . AET. , 

E .. = — —\q\gT eqexp{-E/TAn{-^^) + 1/2 - C] 


d^pd'^x TT 
and the rate for annihilation process is, 


E 


dN 

d^pd'^x 


*^°'^X,\,T'Vexp(-E/T)[M^) - 1 - C] 


TT^ 


4? 


( 22 ) 


Here C=0.577721, is the electric charge of the qnark and the parameter kc is related 

to thermal mass of the qnark in the medinm, 

kl = = 2ml (23) 

with P 

1 A-it 

= + (24) 

In hg.4 and 5, we have shown the photon spectra, as obtained presently, with and withont 
viscons drag at RHIC and LHC energies. It can be seen that photon prodnction is increased 
by a factor of 2 with viscosity turned on. With viscosity turned on, life time of the plasma 
is increased nearly by a factor of 2 and it is natural that photon yield will also be increased. 
However, we have checked that the increase in photon yield beyond pr > 2 GeV is not 
due to increased life time of the plasma, rather, it is due to the slower cooling and lower 
equilibration rate of the of the viscous plasma. The increase in low px yield can be attributed 
to the longer lifetime of the plasma. Also we note that, the photon yield do not depend 
much on the two version of viscosity coefficients used. The yield is marginally increased at 
low pt with viscosity coefficient H, than with viscosity coefficient I, and essentially due to 
the larger lifetime of the plasma. 


B. Lepton pair spectra 


The transverse mass distribution of lepton pairs from nonequilibrium plasma can be 


written as 26 


dN 


dM‘^d?MTdy 


—e^R^ 


f Mt 

J XqXqTdrKoi—) 


T 


(25) 


In £g.5 and 6, we have shown the invariant mass distribution of lepton pairs at RHIC 
and LHC energies. As before, we have shown the results with and without viscous drag. 
Here again, we find that the lepton pair yield in increased by a factor of 2, when the viscosity 
is turned on. As found with the photon spectra, here also, large mass (M > 2 GeV) yield 
of of lepton pairs are essentially due to slower cooling and equilibration rate of the viscous 
plasma, rather than its increased lifetime. The increase in the low mass yield can be said 
to be due to the increased lifetime of the plasma. As with photon spectra, the lepton pair 
yield with two different viscosity coefficients do not differ much. 

The analysis clearly indicate that photon and dilepton yield in QGP are greatly affected 
if dissipative effects like viscosity is taken into account. Large invariant mass lepton pairs 
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and also large pt photons yields are increased by a factor of two, from their ideal connterpart. 
This increase is essentially dne to lower cooling and equilibration rate. The low pt photons 
and the low invariant mass dileptons are also increased, but mainly due to increased life 
time of the viscous plasma. 


IV. SUMMARY 

We have investigated the effect of viscosity on chemical equilibration of parton distri¬ 
bution in relativistic heavy ion collisions. We assumed that after a time Tiso, the partonic 
system achieved kinetic equilibrium. Beyond Tiso expansion (which we assumed to be boost- 
invariant longitudinal only) is governed by the hydrodynamic equations. The approach to 
the chemical equilibration process is then governed by a set of master equations which are 
driven by the two body reactions gg qq and gg <->■ ggg. The partonic matter continue to 
expand and cool due to expansion and chemical equilibration, till the critical temperature 
(Tc= 160 MeV) is reached. Initial condition of the plasma was taken from HIJING calcu¬ 
lations. It was seen that with viscosity turned on, the temperature of the partonic system 
evolve slowly, in comparison to the ideal system. This is expected, as viscosity impedes the 
expansion, making the cooling a slower process. Interestingly, we hnd that the chemical 
equilibration rate is also slowed down in presence of viscous drag. This is in contrast to the 
ideal partonic system behavior, where, lower cooling rate implies a faster equilibration. We 
thus hnd that even though the life time of the plasma is increased by nearly a factor of 2, 
the gluon or quark fugacity could not attain the value that is reached in case of ideal how. 

The ehect of lower cooling and equilibration rate on photon and lepton pair production 
was also studied. It was seen that with viscosity, the photon as well as lepton pair yield 
increased by nearly a factor of 2, both at RHIC and LHC energies. The increase in large 
Pt photons or large invariant mass dileptons are essentially due to the lower cooling rate, 
rather than the large life time of the viscous plasma. The increase in the low pt photon or 
low invariant mass dileptons can be attributed to the to the longer life of the plasma. 

To conclude, viscosity can have a signihcant ehect on parton equilibration process, by 
slowing down the cooling rate and the chemical equilibration process. Its ehect is also 
manifest in the pre-equilibrium photon and dilepton production, raising their yield by nearly 
a factor of 2. 


9 



REFERENCES 


[1] K. Geiger and B. Muller, Nucl.Pliys. B369,600(1992); K. Geiger, Phys. Rep. 
258,376(1995). 

[2] K. J. Eskola and X. N. Wang, Phys. Rev. D49,1284(1994). 

[3] T. S. Biro, E. van Doom, B. Muller, M. H. Thoma and X. N. Wang, Phys. Rev. 
G48,1275(1993). 

[4] K. Geiger and J. I. Kapusta, Phys. Rev. D47,4905(1993). 

[5] K. Geiger, Phys. Rev. D48,4129(1993). 

[6] J. Alam, S. Raha and B. Sinha, Phys. Rev. Lett. 73,1895(1994). 

[7] S. M. H. Wng, Nucl. Phys. A607, 442(1996). 

[8] R. Anishetty, P. Kohler and L. McLerran, Phys. Rev. D22,2793(1980). 

[9] R.G. Hwa and K. Kajantie, Phys. Rev. Lett. 56, 696(1986). 

[10] J. P. Balizt and A. H. Mueller, Nucl. Phys. B289, 847 (1987). 

[11] K. Geiger, Phys. Rev. D46,4965(1992); 46,4986(1992). 

[12] 1. Kawrakow, H. J. Mohring and J. Ranft, Nucl. Phys. A544,471c(1992) 

[13] B. Muller, in Particle Production in Highly Excited Matter, edited by H. Gutbrod and 
J. Rafelski, NATO-ASI Series B( Plenum Press), vl. 303, p.ll. 

[14] E. Shuryak, Phys. Rev. Lett. 68,3270(1992). 

[15] B. Muller and X. N. Wang, Phys. Rev. Lett. 68,2437(1992). 

[16] K. J. Eskola, B. Muller and X. N. Wang, Phys. Lett. B374,20(1996) 

[17] D. K. Srivastava, M. G. Mustafa and B. Muller, Phys. Lett B396,45(1997). 

[18] D. K. Srivastava, M. G. Mustafa and B. Muller, Phys. Rev. G56, 1064(1997). 

[19] P. Danielewicz and M. Gyulassy, Phys. Rev. D31,53(1995). 

[20] A. Hosoya and K. Kajantie, Nucl. Phys. B250,666(1985) 

[21] A. K. Ghaudhuri, Phys. Rev. G51,R2889(1995). 

[22] S. Gavin, Nucl. Phys. A435,826(1985). 

[23] G. Baym, H. Monien, G. J. Pethick and D. G. Ravenhall, Phys. Rev. Lett. 64,1867 
(1990). 

[24] M. H. Thoma. Phys. Lett. B269,144(1991). 

[25] E. Braaten and R. D. Pisarrski, Nucl. Phys. B337, 569 (1990). 

[26] M. Strickland, Phys. Lett. B331,245 (1994). 


10 



TABLES 


Initial conditions characterising the parton plasma at the onset of hydrodynamic evolution 



RHIC 

LHC 

Tjso (fm/c) 

0.31 

0.23 

To(GeV) 

0.57 

0.83 


0.09 

0.14 


0.02 

0.03 
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FIGURES 


FIG. 1. Evolution of temperature, gluon fugacity and quark fugaciy with proper time at RHIC. 


FIG. 2. Evolution of temperature, gluon fugacity and quark fugaciy with proper time at LHC. 


FIG. 3. Transverse momentum distribution of photons at RHIC 


FIG. 4. Transverse momentum distribution of photons at LHC 


FIG. 5. Invariant mass distributions of lepton pairs at RHIC. 


FIG. 6. Invariant mass distributions of lepton pairs at LHC. 
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